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Abstract 

A dielectric composed of a sequence of similar thin layers of alternating refractive index is well 
known to make a good mirror when radiation close to normal incidence has a wavelength within 
the medium of around four times the layer thickness. Such a device — often termed the Bragg 
reflector — is usually introduced to students within the first years of an undergraduate degree, 
and often in isolation from other parts of the course. Here we show that the basic physics of wave 
propagation through a stratified medium can be used to illustrate some more modern developments 
in optics as well as quantum physics; from transfer matrix techniques, to the optical properties 
of cold trapped atoms, optomechanical cooling, and a simple example of a system exhibiting an 
appreciable level of optical non-reciprocity. 
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I. INTRODUCTION 



To paraphrase Griffith^, wave motion can be divided into at least three kinds. These are 
travelhng waves (continuous spectrum); bound states (discrete spectrum); and waves within 
a periodic medium (continuous spectrum with forbidden regions). There are a variety of 
ways in which the third of this number can be introduced. Perhaps the most obvious route is 
through condensed matter physics, where the atomic lattice provides a ubiquitous example 
of a periodic medium (see e.g.'^^. To all intents and purposes such a system has an infinite 
number of periods, where Bloch's theorem approaches exactitude. 

Yet it is instructive to observe the emergence of Bloch's theorem as the number of periods 
of the system is increased. Indeed the motion of waves in a one dimensional periodic, layered 
medium can be solved exactly for any number of layer^^. A neat example where the onset 
of a band gap may be observed can be found in the optics of the Bragg reflectoi^ an 
optical medium composed of layers of two kinds of material with differing refractive indices. 
Possibly the most striking lesson one learns through studying an arbitrary number of layers 
is how quickly the results of Bloch's theorem become applicable as the number of layers is 
increased from unity. 

Not only is the Bragg reflector useful for discussing the emergence of Bloch's theorem, 
but it also unveils some subtleties of propagation in a periodic medium. Dispersion and 
dissipation significantly alter the picture, and dissipation in particular can cause a break- 
down of the distinction between allowed and forbidden bands of frequency. In this work we 
derive the optical properties of the Bragg reflector through an application of the transfer 
matrix formalisnP, and the results are applied to the discussion of the optical properties of 
a lattice of cold trapped atom^, where both dispersion and dissipation are significant at 
the frequencies of interest. As we shall see, the picture becomes quite different when these 
effects are included. 

Further to this we show that the Bragg reflector can be used to introduce the prin- 
ciples of two phenomena of current interest: optomechanical coolin^^ and optical non- 
reciprocitj^^^^. In both cases the high frequency sensitivity close to a band edge significantly 
couples the optical response to the centre of mass motion of the reflector, and thereby pro- 
vides a simple example system where such motional effects can be evident. 
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II. STRATIFIED MEDIA AND TRANSFER MATRICES 



In many cases, knowledge of the transfer matrix reduces the problem of solving a linear 
wave equation in an inhomogeneous medium to that of multiplying together 2x2 matri- 
(3g ji | 8 | i3 | fjj^g wave could be electromagnetic, quantum mechanical, vibrational, fluid me- 
chanical, and so on — so long as the dynamics of the wave-propagation is linear. In the case 
of electromagnetic waves in periodic, layered, planar media the formalism is particularly 
neat. It therefore seems appropriate to introduce the Bragg reflector via the transfer matrix 
formalism. 

The system under consideration is illustrated in figure [T| where we have a dielectric 

(/i = /io) medium composed of homogeneous layers, each with a different value of e = 

eo^^, immersed in background medium of ec = eonl. For arbitrary angles of incidence the 

propagation of waves through the medium depends on the polarization, and in anisotropic or 

moving media, the transfer matrix can become a 4 x 4 object, which significantly complicates 
matterPHm. in 

our case it is sufficient to assume normal incidence onto a medium composed 
of isotropic layers, where the behaviour of the two polarizations is degenerate, and the 
transfer matrix is a 2 x 2 object. 

For a fixed frequency, u the solution to the wave equation within each slice of the medium 
(illustrated in figure [T| is a plane wave, with the spatial dependence of the field being 
proportional to exp {^ik^x) {k^ = n{uj)u)/c is different in each slice). With this in mind, the 
electric field at every point in space is written as a sum of two parts; a right moving part, 
e^~^\x)^ and a left moving part, e^~\x), with the understanding that the full field is the 
sum of the two, E{x) = x^^^'e^^^x) + ''^'•^''(x), where x^^^ are the unit vectors denoting 
the polarization of the waves. The 2x2 matrix, M, that relates these two parts of the 
field across the whole structure is the transfer matrix, and the components of M determine 
the reflection and transmission coefficients of the medium. The only question is how to 
determine M. 

In this case we calculate the transfer matrix associated with getting from the background 
medium, Uc, on the left of figure [ij through the layers, to then again emerge into the 
background medium, Uc, on the right. There are only two kinds of things that happen 
to the light during propagation: it crosses boundaries, and it moves through homogeneous 
media. The full transfer matrix is thus broken up into two kinds of component matrices, an 
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N layers 




FIG. 1: An layer structure is surrounded by a medium with a real refractive index Uc- 
In the rest frame of the medium, we imagine linearly polarized light of frequency, u is 
incident from the left with electric field amplitude, Ei, some of which is reflected as 
(reflectivity, R{u)), and some of which is transmitted as Et (transmittivity, T{uj)). 



^Hnterface^^ matrix, W, that relates the field across a boundary, and a translation^^ matrix, 
X, that relates the components of the field after propagation through a homogeneous region. 
In both cases superscripts are introduced into the notation to indicate the media being 
referred to. The operation of these two matrices is illustrated in figure |2] 

The matrix, W can be found through applying the condition of continuity of the tan- 
gential E and H field^. For normal incidence and non-magnetic media (// ~ /xq), e.g. 
E = EaZ and H = {i/ fiQUj)ydEa/dx, this amounts to the continuity of the electric field 
and its first derivative across the boundary. Applying these conditions on the electric field 



across an interface between index Ua and Ub (see 2a) at x = const., we have 



.{+) 

-'a 



,{-) 
-'a 



,(+) 



-6 ' 



(1) 



and. 



[ei+^ - ei = rib 



(2) 



where we have used the relation, de^^'^ /dx = ii/cajC*^^^ to obtain (^. The above two equations 
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FIG. 2: We consider the total transfer matrix, M, to relate the electric field amplitude of a 
light wave passing from the left to the right of the multilayer stack in figure [TJ In general 

M can be decomposed into a product of 2N + 1 component matrices that are each of two 

kinds denoted by W and X, the actions of which are illustrated in subfigures (a) and (b). 
In (a), the "interface" matrix, W^°'^\ relates the electric field on the immediate left of a 

boundary (here with index n^), to the field on the immediate right (rif,). Meanwhile, in (b) 
the "translation" matrix X^°'\x2 — Xi) relates the electric fields between two points in a 
homogeneous medium (xi and X2), that we here take to have index n^. 



may be cast in matrix form. 





and through inverting the 2x2 matrix on the left, the matrix associated with crossing the 
boundary Ua — ?■ rib is obtained. 



1 / 1 + Ua/Ub 1 - ria/Ub 
1 - Ua/Ub 1 + Ua/rib 





(3) 



The matrix in ([s]) defines the general relation between the electric fields across a boundary 
between two media with different refractive indices. Note that when Ua = Ub, W^""^^ = I2, 
and also that the matrix is transitive, W'-'"'^ ■ W^"''^ = W^''''\ 



5 



The operation of the matrix, X associated with propagation of the field through a dis- 



tance, X2 — X1 within a homogeneous medium is illustrated in figure 2b and is simply obtained 
through multiplying e*^"*"^ and e'-"-' by the phase factor, exp (±ikx(x2 — Xi)), e.g. for a medium 
of index n^. 

The matrices, ([S]) and Q can be used to study light propagation in generic stratified 
media. In the limit of infinitesimally thin slices they can also be used to treat continuously 
inhomogeneous medisP. From here on we consider media with periodically varying layers of 
finite thickness. 



ei^\x2)] 









III. THE BRAGG REFLECTOR: A ID PHOTONIC CRYSTAL 

The transfer matrix approach discussed above is now applied to a stratified medium 
composed of a periodic sequence of identical unit cells (figure |3]), each cell containing two 
isotropic dielectric layers having different complex refractive indices, and thicknesses: such 
a structure is often termed a Bragg reflector. RayleigfP^, and Brag^^ were among the first 
to begin to develop an understanding of wave propagation in periodic media, but it is only 
relatively recently that it has been possible to manufacture three dimensional structures 
with inhomogeneities on the scale of an optical wavelengthpSIZI]_ This, along with theoretical 
work illustrating interesting properties such as inhibited spontaneous emissio d^^*^^ *, photon 
localizatiorP^, and negative refractiorP^, has seen these systems garner the separate general 
designation "photonic crystals". Consequently, in modern terms the Bragg reflector can 
be recognised as an ordered one dimensional photonic crystal, despite the fact that the 
optics of thin films significantly predates this modern worlpSEZI xhe periodic variation in 
the optical response of such media leads to Bragg scattering of light in much the same 



way as for electrons in crystalline solids^ (c.f. sect. Ill A), and this generally leads to the 
appearance of photonic band gaps with a very steep frequency variation of the reflectivity 
and transmissivity profiles around the band-edge. 

The unit-cell of our periodic dielectric multilayer is sketched in figure |3| The relation of 
the electric field at the beginning of cell i to that at the beginning of cell {i + 1) can be 
envisaged as a sequence of W and X operations, as in figure [2] The unit-cell is identified 



cell i-1 



cell i 



cell i+1 




FIG. 3: The unit cell of a Bragg reflector (c.f. figure 1). This is chosen so that it begins 
and ends infinitesimally to the right of the nb — boundary. The relationship between 
the fields in cell i and i + 1 is determined by Ali. The relationship between the field just 
to the right of the rih — )■ ria boundary in one cell, and just to the right of the same 
boundary N unit cells later is instead determined by A^at (cf. (jsj)). For a truly infinite 
periodic structure, the difference in the field across a unit cell is a phase factor e*'*("+''), 
determined by the Bloch wave-vector, k and the crystal periodicity (a + 6). 



as starting and finishing within a medium of index n^, which is equivalent to the following 
sequence of matrices, 

([cos (hb) + iaab sin (hb)] e*''"" —i(3ab sin {kbb)e~'-'^''"- 

iPab sin {kbb)e'''^°'°' [cos {kbb) — iaab sin {kbb)] e'""^"' 

where, = Ua^uj/c, aab = {nl + nl)/{2nanb), and (3ab = {nl - nl)/{2nanb) (recall Ua and 
Ub are complex functions of u). 

The matrix associated with the passage of light through identical unit cells is the A^th 
power of the above unit-cell matrix, AIat = M-i . Given that Ali is a 2 x 2 matrix, we can 
write in a particularly elegant form involving Chebyshev polynomials of the second 
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kincP^. To see how this is done, consider the square of a general 2x2 matrix, A, 

O-li + 012^21 Ctl2(ail + Ct22) 

a2i(aii + 

= Tr[A]A-det[A]l2 (6) 

Having obtained ^ it is clear that it must be possible to write the A^th power of A in the 
form = aNA + bNl2. Multiplying this proposed form on the left by A, and applying 
we obtain the recurrence relations, aN+i = ctArTr[A] + b^, and ^at+i = — aArdet[A], which 
is equivalent to, otv+i = Tr[A]aAr — det[A]aAr_i. Meanwhile, the Chebyshev polynomials 
satisfy a very similar recursion relatioEp^, 



Un+i{z) = 2zUNiz) - Un-i{z) (7) 

this relation holding for both the Chebyshev polynomials of the first, Tn{z), and second, 
Un{z) kind^. The difference between Tjq{z) and Un{z) is in the choice of initial values: 
{To = l,Ti = z}, and {Uq = l,Ui = 2z}. If det[A] = 1, then the coefficients, ajy satisfy 
0, and we find that A^ = Un-i{z)A - Un-2{z)12, with z = Tr[A]/2. We must choose 
the Chebyshev polynomial of the second kind because we must obtain (|6| when N = 2. We 
assume det[AlAr] = 1, which for our purposes is always the cas d^^*^^^ (see (34) below). 
The A^-cell matrix thus takes the form, 



where, 



MN = UN^liz)Mi-UN-2iz)t2, (8) 



z = -Tr (A^i) = cos {kfyb) cos {kaa) — a^a sin {ki,b) sin (fcatt). (9) 



with All given by ([s]), and where the Chebyshev polynomials may be computed from their 
relation to the trigonometric functiond^ Un{z) = sin [(A^ + 1) arccos(z)]/-\/l — z'^. 

The transfer matrix, M, for the situation illustrated in figure [T] is A4.N sandwiched 
between two matrices associated with light passing from the background index, Uc, into the 
multilayer. 



M 



^ niii mi2 

ym21 17122, 



W^"'^ ■Mn-W^^''\ (10) 
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where we apply W^'^^^ •W^"-'''' = W^°-'^\ The relationship between the electric field amplitudes 
on the left and right of figure [l] is finally given by Cr = M • Cl, where eL,R here refers to the 
two row vector composed of the left and right going parts of the field. 

Refiection from, and transmission through the system can now be calculated: light of 
amplitude Eq is incident from the left, and the transfer matrix equation reduces to, 

mil 

17121 

from which follows expressions for the refiectivity and and transmissivity. 





i?4 



det[M] 



m22 



m2i 



m22 



(11) 



assuming that ric is real. Similarly, for incidence from the right we have, 

mil mi2 
17121 m22 





which gives. 



1 



1^22 

mi2 



m22 



(12) 



Consequently (11) and (12) are not equal in general. This is not too surprising. It is 



straightforward to imagine multilayer systems where identical refiectivity for incidence from 
the left and the right would not be expected. One example would be a highly refiective 
surface backed by a thick layer of strongly absorbing material: incidence onto the front 
surface would yield a refiectivity close to unity, while incidence onto the reverse would yield 
a much lower value. The situation is not so clear in the case of the transmission coefficient, 
as one requires det[7Vf] 7^ 1 in order for ^ T_: such a system would be optically non- 
reciprocaP^^, which is a non-trivial property in electromagnetism. We shall return to this 
issue in section IIV CI 

For our specific system the transfer matrix is (10). Inserting ([s]), ([S]), and ^ into (10) 
gives the explicit expressions for the refiectivity and transmissivity of this system. The 



elements of the transfer matrix are found to be 



where, 



mil = Un^i{z) [z + iOacC - iPacPab COS [katt) sin {kbb)] - Un-2{z) 
mi2 = -Un-i{z) [Pab sin {kbb) (sin (/c^a) + iaac cos {kaa)) - i(3acC] 
17121 = -Un-i{z) [f3ab sin {kbb) (sin {kaO) - iaac cos {kaa)) + if3acC] 
"^22 = Un^i{z) [z - iaacC + 'iPacPab COS {katt) sin {kbb)] - Un-2{z) 

( = sin {katt) cos {kbb) + aba cos {kaa) sin (/i;fe6), 



(13) 



(14) 



The results (13) can then be combined to give the desired reflection and transmission coef- 



ficients. The transmissivity is found to be equal for incidence from the left and right, 

T 2 



Un-i{z) [z - iaacC + iPacPab COS {katt) slu {kbb)] - Un-2{z) 
Meanwhile the reflectivity is not, 

Un-i{z) [f3ab sin {kbb) (sin {kaa) - iaac cos {kaa)) + i/3acC] ^ 



(15) 



Un-i{z) [z - iaacC + iPacPab COS (fc^a) sin {kbb)] - Un-2{z) 
Un-i{z) [f3ab sin {kbb) (sin {kaa) + iaac cos {kaa)) - i/3acC] 



(16) 



Un^i{z) [z - iaacC + iPacPab COS (fc^a) sin {kbb)] - Un-2{z) 

unless the medium is lossless, in which case it is evident that mi2 = "^21. The reason for this 
difference is that, when incident from the left, light enters the medium initially into a region 
of index equal to n^, and initially into a region of index equal to Ub when incident from the 
right. For thin layers with small dissipation, the difference is negligible. The expressions for 
i?+ and T+ are plotted in figure 111 for a multilayer of 10 unit cells. 



A. An optical Bloch theorem 

Any real multilayer system will have a finite number of layers; a Bragg reflector may 
only have tens or hundreds. It is therefore perhaps surprising that we can glean some of the 
behaviour illustrated in the previous section — culminating in figure |4] — from the case of an 
infinite number of layers. Nevertheless, Bloch's theorem can be used to accurately predict 
the region of high reflectivity in the central region of figure |4] (20 layers). The reason for its 
applicability is that for such a range of frequencies, the wave is rapidly extinguished from 
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FIG. 4: Transmissivity Tj^iiS) and reflectivity -R+(u;) as a function of frequency for a 
transparent Bragg reflector (pliotonic crystal). The system is composed of 10 unit cells 
where all layers have a real refractive index, na = 1.45 and rih = 2.1. The thicknesses of the 

layers are a = 2.60 x 10^^ m and b = 1.76 x 10^'' m. The spectra contain a stop-band 
region where propagation is forbidden, corresponding to the plateau in the central region 
of or to the grey shaded region where \z\ > 1 (inset). The amplitude of the oscillations 
increases with the index contrast from = 1 — )■ = 1.45 at the edge of the crystal, but 
this does not affect the stop band structure. Such an increase in oscillation amplitude can 



be important (see section IV). Inset: For an infinite array of unit cells no mode can 
propagate in the frequency range shaded in grey {\z\ > 1, corresponding to the photonic 

band gap. 



the Bragg reflector so that the field amplitude is very close to zero at the exit interface, and 
we can consider the number of unit cells to be infinite with little error. 

In an optical system with translational symmetry along an axis, x, the eigenstates, E 



11 



transform as E{x + S) = E[x)e^^^ , due to the fact that the translation operator commutes 
with the operators within the electromagnetic wave equation. Bloch's theorem is the ana- 
logue of this relation for the case when we do not have continuous translational symmetry, 
but instead an infinite periodic structure, with discrete translational symmetry. In this case, 
if the spatial periodicity is a + 6, then the eigenmodes obey, 

E{x + a + h) = E{x)e'^^^+'''^ (17) 



The quantity, k is known as the Bloch wave-vector (see figure |3] for illustration). In general 
K is complex, and designates a normal mode for light waves that propagate by coupling to 
the periodic medium. The Bloch wave-vector is the analogue of the free space wave-vector 
for the case of a discrete translational symmetry, and tells us the large scale variation of the 
field as one moves across one or more unit cells. In terms of the transfer matrix formalism 
outlined in the previous two sections, this statement is, 

Mi-e = e'"('^+'')e, (18) 



Equation (18) can be seen as analogous to the electronic Bloch theorem, where the unit 
cell transfer matrix, Ali is playing the role of the translation operator, and the electric 
field e is cast as the wave-function. Finding the relation between the frequency, u, and the 
Bloch- vector, k amounts to finding the eigenvalues, A = e*'^^""'"^) of the unit-cell matrix, Adi. 



These eigenvalues can be found in the usual manner, taking the right hand side of (18) to 
the left and setting the determinant of the resulting matrix to zero. 



det {Ml - AI2) 



Mu - A M12 
M21 M22 - A 



0, 



or. 



A^ - A {Mil + M22) + {M11M22 - M12M21) = A^ - 2zA + 1 = 0, (19) 



where we have used the fact that det[A4i] = 1. The quadratic equation on the right of (19) 
has the two solutions, 

A = ^ ± = e^'^('^+'') (20) 

These solutions determine the dispersion of the Bloch modes, i.e., the function k,{u). It can 
be seen directly from ^ that z = Tr[Ali]/2 is real for media with real refractive indices 
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and Ub- In this case the criterion for allowed propagation through the multilayer becomes 
quite simple, for when 2; < 1 we have, 

which can be satisfied with a real Bloch wave-vector: k = ± arccos(2;)/(a + b). Meanwhile, 



when z > 1, the right hand side of (20) is real, with one root greater, and one less than 
unity. This leads to a complex value of k. Such solutions have divergent behaviour at 
infinity, and are therefore not allowed modes of the truly infinite system. The frequency 
intervals where z > 1 correspond to the forbidden gaps where the real part of n is restricted 
to be either (for gaps opening at the centre of the Brillouin zone) or ±7c/{a + b) (opening at 
the Brillouin zone boundaries). While the former solutions with real k describe light waves 
that propagate in the medium, the latter ones with complex k, can be interpreted as surface 
excitations {Tamm states). In a bounded medium these are evanescent waves that undergo 
extinction inside the medium in the direction normal to the boundary. The decay is due to 
strong Bragg reflections, a reversible process that does not lead to a loss of energy from the 
field, making reflection within the gap strictly unity. 

In a dissipative periodic medium the situation becomes rather intricate: z is complex. To 
see how this modifies the situation, consider the case where a small degree of dissipation is 
present: z = cos{ip) +iT], where both (f and rj are real, and where rj <^1. To first order in 77, 

6'''^''+^'^ ~ (1 ± -^^] e^'^ (21) 



The presence of dissipation shifts the modulus of (21) away from unity. Therefore k is in 
general complex when dissipation is present, even for modes falling within an allowed band. 
This has the simple interpretation that in an absorbing periodic medium, the propagation 
of Bloch waves is damped. This is a dissipative process, which prevents reflection from 



being unity. It is also clear that the forbidden gaps are significantly modified (see figure 6a) 



When the real part of k is close to or 7r/(a + b) — where the forbidden band would begin in 



the absence of dissipation — then sm{(p) ~ 0, and the right hand side of (21) has a modulus 



significantly different from unity. Hence the distinction between allowed and forbidden bands 



becomes blurred (see section IV A). For large degrees of dissipation the band-gap may even 
disappear. 
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IV. SOME EXAMPLE APPLICATIONS 



So far the discussion may have been rather famihar. Although we have given a treatment 
of the Bragg reflector in terms of transfer matrices, which might not be the standard presen- 
tation for undergraduates, the formulae are mostly well-known. We shall now demonstrate 
how we may use these results to discuss some examples that are relevant to contemporary 
experimental set ups. 

A. A Bragg reflector of cold atoms loaded into an optical lattice 




FIG. 5: Neutral atoms can be trapped within the electric field of two counter-propagating 
red-detuned laser beams. The electric field induces an atomic dipole moment, which 
interacts with the field, pushing the atomic centre of mass towards high field intensity. 
When sufficiently cooled these atoms can be trapped, and form a lattice. A weak probe 
beam interacts with the optical lattice in much the same way as the Bragg reflector 

illustrated in figure j3j 

A neutral atom subject to an external electric field will develop an electric dipole mo- 
ment, d. The atom will then be subject to a force, F = dp/dt, due to interaction of the 
electromagnetic field with this induced polarization, and this force can be used to trapa 
cloud of such atoms within a confined region of spac^^. In the situation illustrated in fig- 
ure |5| the interference pattern created by two counter-propagating laser fields creates an 
array of many thousand of such traps, often termed an optical lattic^. Once cold atoms are 
loaded into such a system then one has an array of interacting sites, each containing some 
average number of atoms per site, where the interaction strength of the atoms can be tuned 
via the laser field intensity. Within such a system one may create analogues of many body 
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phenomena in condensed matter physic^. 

Consider a monochromatic field in which the induced atomic dipole moment can be 
written to first order in the electric field as, d = a{u))E{r), where a{uj) is the single-atom 
polarizability (assumed real and positive, as appropriate for a red detuned optical lattice), E 
is the electric field amplitude (in this case a three dimensional vector), and r is the position 
of the atomic centre of mass. The force on the centre of mass of each atom is, to first order 
in the size of the system. 



dp 



j [pE + j X B]dV ^ j xp{x + r)dV ■ VrE{r) + J ji^ - ^)dV X B{r) 



where the integral is over the volume of the system of charges. We assume that the velocity of 
the centre of mass is slow enough that we can neglect contributions to the force of order f/c, 
and we also neglect the magnetic structure of the system. Identifying, d = j xp{x + r)dV , 
and d = J j{x — r)dV^, the force on the centre of mass can be written as, dp/dt = 
{d • ^r)E{r) + d X B. For the particular case of an induced dipole moment this then 
becomes, 

^ = la{u:)V\E{r)\' + a(a;)| {E x B) (22) 

Figure [5] illustrates two counter-propagating optical beams that create the electric field that 
the atom encounters. Averaging the motion over a timescale significantly longer than the 



period of the optical field, the final term on the right (22) becomes negligible. The average 
force is then such that the atomic centre of mass will be pulled towards regions of high electric 
field intensity. In figurejsjthe red line indicates the effective potential, V^g = —^a{u)\E{r)\^ ^ 
in the minima of which we can imagine atoms collecting once their centre of mass motion 
has been cooled significantly by other mean^. 

Not only is the optical lattice useful as a model system, but it also has interesting optical 
properties of its owEpESl \Ye can treat this system using the formalism developed in the 
previous sections (see figure [s]). In the case of a one dimensional optical lattice, where the 
atomic motion is only weakly constrained in the y — z plane, the system is essentially a 
slightly unusual Bragg refiector (an 'atomic Bragg mirror'). The refiector is approximately 
composed of layers of atomic clouds (width a, index n^), and layers of vacuum (width 6, 
riij^ric ~ 1). 

Suppose that a weak beam of light probes the atoms in the lattice; what are the reflection 
and transmission coefficients? When the frequency of this probe beam is close to an atomic 
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W(Hz) 



(a) 



(b) 



FIG. 6: Transmissivity T+(a;) and reflectivity R+{uj) as a function of frequency for a 
dissipative Bragg reflector, (a) Tlie dissipative analogue of the plot shown in figure |4j All 
parameters are identical, but with the addition of an imaginary part to the refractive 
index, nf,. As Im[nb] is increased, the contrast of the fringes around the stop band is 
decreased, and the shape of the reflectivity profile in the stop band changes, (b) The array 
of unit cells now describes the index profile of cold atoms trapped within an optical lattice 
(A/" = 5.7 X 10~^). The probe beam is assumed to be close to the frequency of the 



Si/2 -P3/2 transition of '^^Rb, where the index can be approximated by (23), with 
uo = 2.412 X 10^^ Hz, 7e = 3.7699 x 10'^ Hz, and we use the notation, 6 = {uq — oj)l^e- The 
trap has = 5.4 x 10^ periods, and the periodicity is such that a = 1.94807 x lO^'^m, and 
b = 3.70873 X 10^^ m. The blue and dashed lines illustrate the effect of changing the line 

shape parameter, (. 



resonance, the atom may be approximated as a two level system. The index of such a cloud 
of two-level atoms can then be calculated using the Maxwell-Bloch equationd^, taking the 
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approximate form, 



Ua (w) 



'1 + 



(Wo - Uj)/le - <' 



(23) 



In (23) 7e is the linewidth of the excited state of the trapped atoms, is a parameter which 
modifies the shape of the resonance (C = 1 being the theoretical result emerging from the 
two-level approximation), J\f is the scaled atomic density, A/{kQV) {A/V atoms per unit 
volume), and ko = uo/cis the resonant wave-vector associated with the frequency of atomic 
resonance, uq. 
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FIG. 7: (a) The absolute value of the right hand side of (20). The deviation of this 
quantity away from unity indicates the magnitude of the imaginary part of the Bloch 
wave-vector, (b) The real and imaginary (inset) parts of the refractive index computed 



from ( 23 ) . In both cases the parameters and variables are as in figure 6b 



Figure 6b shows the reflection and transmission coefficients calculated from (15-16) 



for the case of trapped ^''Rb. We can see that close to the resonant frequency (5 = 0), 
where the dissipation and dispersion are both most pronounced, the optical response is 
quite complicated, with rapid oscillations evident on the scale of the 7e. Either side of 



6 = two bands of high reflectivity appear. Figure 7a is a plot of the absolute value of the 



right hand of (20), illustrating the origin of these two bands. For negative detuning away 



from resonance — where flgure 6b shows nearly total reflection — the Bloch wave-vector has 



a signiflcant imaginary part even though the absorption is small (inset of flgure 7b ) . There 
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is also a stop band evident for positive detuning at around S ~ 250. This is not completely 



evident in figure [6b} where the reflectivity peaks at 0.6 for 5 ~ 150, but does become 
conspicuous for larger values of A^. 

This example system illustrates that when composed of resonant elements, the properties 
of the Bragg reflector can be quite subtle, and not at all like the textbook examples. Fur- 
thermore it shows that even when the layers have an index contrast of less than a percent 
away from unity, it is still possible to create an almost perfect reflector. This is so long as 
one can construct a system with N large enough, which is quite possible with an optical 
lattic^^SES] xhese atomic Bragg mirrors can be devised to re-direclP^, confind^, or amplifj^ 
radiation. 



B. Optomechanical cooling and heating 



X 




FIG. 8: A Bragg reflector is allowed to move along the x axis, confined by a pair of springs 
attached to fixed columns, and performing an oscillatory motion around the point x = xo- 
When equal amplitude, monochromatic light is incident onto both sides of the reflector, 
radiation pressure can be used to add or remove kinetic energy from the motion. 
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The electromagnetic field can add or remove energy from a mechanical degree of free- 
dom. In much the same way that atomic motion can be cooled within an optical lattice, 
macroscopic mechanical motion can be reduced to the point where it enters regime where a 
quantum mechanical description must be usecP^'^. For example, in the case of a mechanical 
oscillator this would be when the total energy in the centre of mass motion becomes com- 
parable with the spacing of the energy levels. Recent experimental work has succeeded in 
reducing the vibrational energy of macroscopic degrees of freedom to the quantum mechan- 
ical ground stat^. We consider the system illustrated in figure [s] to demonstrate the basic 
principles (if not the detailed experimental set-up) involved in optomechanical cooling. This 
is inspired by the fact that a similar system has been recently proposed for this purpos#^'^. 

We imagine monochromatic radiation is incident onto both sides of a Bragg reflector, and 
due to the reflection from, and absorption into the body, it will exert a force. We calculate 
the optical force on the reflector by firstly examining the average balance of energy and 
momentum for the case of a single photon impinging on the structure at rest. In the rest 



frame of the reflector, we can use the reflection and transmission coefficients (15-16), where 
they are properly defined. For a photon incident from the left with momentum hk+x and 
energy huj+, the difference in optical energy and momentum before and after interaction 
with the reflector would be found to be, on average, 

AE+ = hiu+ [T+{cu+) + R+iuj+) - 1] 

Ap+ = hk+ [T+(cj+) - R+{cu+) -l]x (24) 
Similarly, for incidence from the right with momentum —hk-X and energy hw-, 

AE^ = hcu^ [T-{u-) + R-{u-) - 1] 

Ap_ = -hk_ [T_(a;_) - i?_(a;_) -l]x (25) 

Assuming V± photons per second are incident from the left and right respectively, we can 
equate the rate of energy and momentum being lost from the fleld to that taken up by the 
medium. The rest frame force on the reflector is then, 

^ = h{V+k+ [1 + R{uj+) - T{uj+)] - V-k_ [1 + i?(a;_) - T{uj^)]}x (26) 

and the energy transfer (heating) is, 
dE' 

^ = h {V^u^ [1 - R{uj^) - T(a;+)] + V^uj_ [1 - R{uj_) - T(c._)]} (27) 
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assuming i?+ ~ _R_ = R and T+ = T_ = T. 

Equal amplitude monochromatic radiation of frequency u incident onto both sides of the 
reflector is assumed in the laboratory frame. To first order in the velocity, the rest frame 
frequency and wave- vector are given in terms of laboratory quantities by, u)± = w(l=Fx/c) = 
ck±. Furthermore, in the laboratory, the number of photons sent towards the medium per 
second, Vq is equal on both sides of the reflector, and this also transforms as a frequency, 
V± = Vo{l T i / c) . Expressions (26-27) are then. 



dt 







m 









dR{uj) dT{uj) 



2 [1 + R{lo) - T{lo)] )■ X 



(28) 



and. 



dt 



2hujVo [1 - Riuj) - Tito)] 



(29) 



where the reflection and transmission coefficients have been expanded as, e.g. R{u)±) 



R{ui) =F {xui/c)dR{uj)/duj. Finally, transforming (28-29) from the rest frame into the labo- 
ratory frame, using dp^^i/dt ~ dp'-^^Jdt + {x/c^)dE'^-^/dt, the optical force on the reflector 
due to the field is then of the form, 



dp^ 



dt 



1 + 3R{uj) - T{uj) + u 



dR{u) _ dT{uj) 
du du 



X 



(30) 



Expression (|30|) is the optical contribution to the force on the medium, as observed in the 

dR{u) dT{u) 



laboratory frame. Defining the following quantity 

huVo 



1 + 3R{u) - T{u) + u 



du 



du 



(31) 



the equation of motion for the centre of mass of the Bragg reflector in the presence of the 
optical field becomes, 

(32) 



m—— + 21 — — h mil (x — Xq) 
dt^ dt 



where m is the mass of the refiector, and Q is the frequency associated with the restoring 



force of the springs (see figure [8j). The equation of motion, (32) is of the form of a damped 
(r > 0), or amplified (F < 0) harmonic oscillator. The damping coefficient is proportional 
to the incident optical power, and depends on the refiection and transmission coefficients 



and their derivatives with respect to frequency. The solutions to (32) have the form, 

xit) =Xo + Aoe-r*/'"e±*^Vi-rV™^f^^t 
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FIG. 9: The damping coefficient, F in units of huVo/c^, for a Bragg reflector in a 
frequency regime around the edge of a stop band. The parameters are as in figure |4| but 
for the case of = 50 layers, where the stop band is more defined. 



where Aq is the distance of the centre of mass from xq at t = 0. If the medium is only 
weakly dispersive, then the derivatives of the reflection and transmission coefficients within 



(31) may be neglected: then F is strictly positive. Without such terms it is evident that the 
rate of damping is ordinarily very small. For instance, an optical field, u ~ 10^^ Hz, and 
a reflector of a very small mass^-^', m ~ 10~^^kg, yields, F/m ~ lO^'^^Vo- For a moderate 
damping rate of F/m ~ 10~^ we would thus require Vq ~ 10^*^, corresponding to a laser 
power per unit area of lOW/^, where A is the cross sectional area of the reflector. In such 
a field the heating due to absorption alone would be enough to melt the medium. 

Conversely for a Bragg reflector, dispersion is not negligible. Figure |4] shows there are 
strong oscillations in the reflection and transmission coefficients, contributing significantly to 



F via the frequency derivatives in (31 ). The motion can then be either amplified or damped, 
depending on the sign of the derivatives of the reflection and transmission coefficients. For 
small dissipation, and increasing reflectivity with frequency we have damping (F > 0), while 
for decreasing reflectivity with frequency, we have amplification (F < 0). Figure [9] illustrates 
c^r /hwVo, demonstrating four orders of magnitude increase in the absolute value of the 
damping coefficient. The required optical power density is correspondingly four orders of 
magnitude smaller. The above effect is the basic principle of optomechanical cooling. 
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It is worth mentioning that while the Bragg reflector is highly dispersive, a typical exper- 
iment in optomechanics involves an object that has relatively weak dispersion. However, the 
object sits within a high quality factor optical cavity. The frequency of incident radiation is 
close to a cavity resonance, which is the source of the strong dispersion in the equivalent of 

C. Optical non-reciprocity 

Reciprocity — the fact that the same readings will be obtained when the positions of a 
source and detector of waves are interchanged — is a rather general property of wave propa- 
gation. For the purposes of performing signal processing, or computation with photons, one 
might like to realise an optical diode, allowing total light transmission in the forward direc- 
tion, and inhibiting (over some bandwidth) propagation in the backward direction. Partly 
for this reason, current efforts are going into the realisation of non-reciprocal optical devices. 




FIG. 10: In the laboratory frame a Bragg reflector moves with velocity v along x. Light of 
frequency u, incident on the multilayer structure, is seen to exhibit signiflcantly different 
reflection and transmission characteristics, depending on whether it impinges from the 

right {-R_,T__}, or from the left {i?+,T+}. 

Non-reciprocal devices are often composed of materials exhibiting magneto-optical ef- 
fect^l^, but these are difficult to integrate into an optical circuit. Nonlinear optics can 
also be usecP^'^, but good performance is then restricted to a specific power range. These 
two methods are by no means the only options, and for reference we note that several other 
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mechanisms for breaking reciprocity have been explorecP^"^ some of which have been re- 
ahzed experimentallj^. Here we illustrate the pervasiveness of electromagnetic reciprocity 
within the transfer matrix formalism of section |Tl} A moving Bragg reflector provides a 
simple example of a system exhibiting significant non-reciprocity (due to the motion, which 
breaks time symmetry). Indeed close to the stop band edge its properties can be much like 
that of an optical diode. 

Consider propagation through an arbitrary layered medium, as described in section |TT[ 
and illustrated in figure [TJ If the medium is reciprocal, then the transmission from the 
background index on the left of the medium to that on the right should be the same as 
transmission from right to left. We shall define a simple measure of the non-reciprocity of 
our layered medium to be, 

AT = T+{uj) -T4uj). (33) 

where the transmission coefficients are those measured in the laboratory. This measure 
obviously doesn't capture everything that the full definition of non-reciprocity does (see 
e.g.'^, but should be in qualitative agreement. Comparing (11) and (12) it is evident 



that reciprocity as defined by (33) requires det [iVf] = 1. For normal incidence onto a non- 
magnetic medium, an arbitrary transfer matrix is equal to a (possibly infinite) product of X 
and W matrices given by (|3]-[4]), and the determinant of such a product equals the product 
of the determinants, 

N N 

det[M] = JJdet[X«(rfi)] JJdet[iy(^'^'+^)] 

i=l j=0 
N 

= TT^ = — (34) 



where we have assumed an N layered medium, each layer having index rij and thickness di 



As defined here, (33) is applicable only to the case where uq = n^+i, and thus T+ = T_. One 
may well ask about the case where we examine the electric field within the medium rather 
than outside, where the above analysis shows that det [iVf] ^ 1. Initially it appears that this 
case violates reciprocity. However this conclusion is incorrect, for the fundamental measure 
of reciprocity is about the field produced by a current of fixed amplitude, as registered by 
a detector. For cases where the source and detector are in different media, we must modify 



the simple measure (33). 
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FIG. 11: Non reciprocal response, (33) for an atomic Bragg mirror set in motion. The 



velocity is 1ms ^ (solid red) and 5ms ^ (black dashed). The parameters are those given in 

figure 



6b 



For the case illustrated above we have shown reciprocity to be necessarily embedded 
within the 2x2 transfer matrix formalism. To break reciprocity we must do something 
to break this formalism. Perhaps the simplest way to do this is to take a medium that 
has been set into motion, as in figure [lOj We cannot then simply imagine that uj is the 
same throughout the system. Refiection causes the frequency of the waves to become u:± = 
uj{l =F f/c), and this will depend on the direction of incidence. Therefore, incident left and 
right going waves of frequency u will generate refiected and transmitted waves of frequencies, 
u± and u: for two input channels we have four output channels. The 2x2 formalism is 
therefore not applicable in the frame where the medium is in motion. 



To calculate (33) for our system, we first examine propagation in the rest frame. For 



E'(x') 



incidence from the left, the rest frame field on either side of the medium is 

E+z [e^^+^'Z^ + r+(w+)e-*'^+^''/'=] x' < 
E+i;t(u;+)e^'^+^'/^ x' > L 

and for incidence from the right, 

E_zt(w_)e-^'^-^'/^ x' < 

E_z \e-^^'^'l^ + r_(u;_)e''^-^'/^] x' > L 



(35) 



E\x') 



(36) 
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where the medium is taken to be of length L. Expressions (35 36) are transformed into the 



laboratory frame using the first order transformation, E = E' — v X B', 



E(x) 



Eoz [e^"^^/^ + (1 - 2f /c)r+(w+)e-''^^/^] x < vt 



iujx/c 



and. 



E(x) 



Eozt{u+)e 



Eoz{l + v/c)t{uj^)e 



(37) 



-iujx/c 



X > L + vt 



X < vt 



Eoz{l + v/c) [e"*^^/" + (1 + 2t;/c)r_(w_)e^'^^/^] x > L + vt 



(38) 



where E^ = Eq{1 — v/c) and E_ = Eq{1 + v / c), the notation being as in section IV B The 
difference in laboratory frame transmissivities for left and right incidence can thus be inferred 



from (37-38) and simply involves the difference in the rest frame transmission coefficients 



for two different Doppler shifted frequencies, 



AT = T. - T_ 



\t{u{l-v/c))\^ -\t{u{l + v/c))f 



(39) 



To first order in the velocity of the medium, (39) is 



AT 



2ujv dT{uj) 
c du 



(40) 



where T{uj) is the rest frame transmissivity. If dispersion is moderate then this measure 
of non-reciprocity is negligibly small for realistic velocities (say T changes 0.1 over 10^^ Hz 
then for optical frequencies and v ~ lms~^, AT ~ 10^^). However, when the rest frame 



transmission coefficient depends strongly on frequency, then (39) may be large: again the 



Bragg reflector illustrates a prime example of this. Figure 11 gives a plot of AT for an atomic 
Bragg mirror (section IV A) set into motiorJ^ at lms~^ and 5ms~^. Due to the strongly 



dispersive optical response of the lattice (figure 6b) the non-reciprocity of this system as 



determined by AT can be in the tens of percent for velocities of metres per second. For 
the same velocity regime, an effective optical diode exhibiting nearly 100% isolation would 
require only an order of magnitude increase in the dispersion. This could be realised in a 
three level systems driven by an optical field, as used to achieve electromagnetically induced 
transparency^^*^. 
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V. CONCLUSIONS 



We have given an overview of the properties of the Bragg reflector using a simple appli- 
cation of the transfer matrix formalism. The basic physics of this system was found to be 
convenient for discussing several topical areas of optical physics. We used it, in particular, 
as the basis for a discussion of optomechanical cooling and optical nonreciprocity effects 
which can all be observed by using Bragg mirrors made of ultra-cold atoms. 
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noise within the optical field. This extra fluctuation serves to counteract the cooling process. 
An atomic Bragg mirror may be set into motion by relatively detuning the two laser beams in 
figure [5] so that the interference pattern moves along at a velocity determined by the detuning. 
Alternatively the lattice may be given a sudden jolt so that the atoms slosh about within the 
wells. In this latter case the velocity of the medium is a periodic function of time. 
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